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FREE-WAVE  PROPAGATION  RELATIONSHIPS  OF 
SECOND-ORDER  AND  FOURTH-ORDER  PERIODIC  SYSTEMS 


1.  INTRODUCTION 


Periodic  systems  are  studied  and  modeled  for  use  in  different  structures  for  various 
rationales.  Some  examples  are  reinforeed  ship  hulls,  railroad  tracks,  aircraft  fuselages,  and 
trolley  lines.  Analysis  of  these  systems  has  been  ongoing  for  many  years.  Some  of  the  analysis 
ean  be  extremely  eomplex  (reference  1)  and  can  produee  a  solution  to  a  speeifie  problem,  and 
some  of  the  analysis  ean  be  simple  (reference  2)  and  produces  insight  into  system  behavior. 

These  problems  are  sometimes  formulated  to  determine  the  free-wave  propagation  of  a  structure 
and  sometimes  to  calculate  the  forced  wave  response  to  a  speeifie  type  of  load.  For  simple 
mechanical  systems,  these  problems  can  be  broken  down  into  two  speeifie  classes: 

1 .  Seeond-order  systems  (such  as  bars  and  strings)  where  the  governing  differential 
equation  contains  a  seeond-order  derivative  with  respeet  to  the  spatial  variable. 

2.  Fourth-order  systems  (sueh  as  plates  and  beams)  where  the  governing  differential 
equation  contains  a  fourth-order  derivative  with  respeet  to  the  spatial  variable. 

Models  of  seeond  order  systems  begin  with  the  wave  equation  in  spaee  and  time.  These 
systems  can  beeome  periodic  with  the  addition  of  stiffeners  (springs)  and/or  masses  (beads). 
These  systems  have  been  studied  with  multiple  sets  of  stiffeners  and  masses  (reference  3)  to 
determine  the  dynamie  response.  The  tensioned  string  with  periodic  stiffeners  has  been  analyzed 
in  the  frequency  domain  with  moving  harmonic  forces  (reference  4)  and  a  suddenly-applied 
eoneentrated  foree  (reference  5).  These  papers  discuss  the  stability  of  the  system,  particularly 
with  respeet  to  the  value  of  the  stiffeners,  although  there  is  no  extension  of  the  problem  into  the 
wavenumber  domain.  The  response  in  a  nearly-periodie  beaded  string  has  been  analyzed  using  a 
matrix  method  (reference  6).  A  model  and  experiment  of  a  nearly-periodie  pendulum  system  has 
also  been  formulated  (reference  7)  with  an  emphasis  on  Anderson  localization. 

Models  of  fourth-order  systems  typically  begin  with  the  Euler-Bemoulli  beam  (or  plate) 
equation  in  space  and  time.  These  systems  ean  beeome  periodic  with  the  addition  of  stiffeners 
(springs)  and/or  masses  (beads).  In  reference  8,  a  fourth-order  system  was  studied  for  free-wave 
propagation  and  a  distributed  pressure  load;  a  spatial  domain  solution  was  constructed  and  a 
contour  integral  was  evaluated  to  transfer  the  expression  into  the  wavenumber  domain.  This 
work  contains  an  expression  for  free-wave  propagation  in  the  reinforeed  plate.  The  work  in 
reference  8  was  extended  to  stiffened  plates  with  line  and  point  forces  (reference  9).  Work  has 
also  been  completed  on  the  modeling  of  plates  that  have  sectional  aperiodicity  (reference  10). 

The  present  report  develops  a  relationship  between  the  determinant  of  a  periodic  system 
dynamie  matrix  and  an  analytical  expression.  Speeifieally,  the  ease  of  free-wave  motion  is 
investigated  in  the  wavenumber  frequency  domain.  Two  separate  methods  are  used  to  solve  the 


periodic  problem.  First,  a  series  solution  to  the  problem  is  generated  and  this  is  converted  into  a 
closed  form  analytical  expression  using  infinite  series  of  residues.  Next,  a  matrix  solution  is 
formulated.  By  equating  functions  of  the  two  methods  for  free-wave  propagation,  an  analytical 
expression  is  derived  when  the  determinant  of  the  system  matrix  is  zero.  This  expression  is 
useful  for  understanding  the  behavior  of  these  diagonally-indexed,  full  matrices.  This  method  is 
applied  to  both  a  second-order  and  fourth  order  system.  The  origin  of  each  individual  problem  is 
discussed,  and  two  numerical  examples  are  illustrated. 


2.  SECOND-ORDER  SYSTEM 


The  motion  of  a  reinforced  string  is  a  second-order  system  that  is  governed  by  the  wave 
equation.  If  the  stiffeners  reinforcing  the  string  are  equally  spaced,  and  if  the  string  has  infinite 
spatial  extent,  then  the  partial  differential  equation  modeling  this  motion  in  the  spatial-time 
domain  is 
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where 


•  u(x,/)  is  the  transverse  displacement  in  the  ^-direction, 

•  x  is  the  spatial  location  on  the  string, 

•  /  is  time, 

•  T  is  the  tension, 

•  p  is  the  mass  per  unit  length, 

•  f(xj)  is  the  external  load  on  the  string  per  unit  length, 

•  K  is  the  stiffness  of  each  reinforcing  member, 

•  L  is  the  distance  between  adjacent  stiffeners,  and 

•  5  is  the  Dirac  delta  function. 


Using  a  point  load  that  is  harmonic  in  time,  and  transferring  equation  (1)  into  the 
wavenumber-frequency  domain  yields  (reference  3) 
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where  \v(k)  is  the  transverse  displacement  in  the  >’-direction,  k  is  the  wavenumber,  Fn  is  the 
magnitude  of  an  applied  point  force  atx  =  0,  and  k,  is  the  free  wavenumber  of  the  unreinforced 
string  and  is  given  by 
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where  co  is  the  frequency  and  c,  is  the  free-wave  speed  of  the  unreinforced  string.  Note  that  the 
free  wavenumber  implicitly  contains  the  frequency.  The  solution  to  this  equation  has  been 
previously  solved  and  is  given  by 
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where  k„  =k  +  2/r nl  L  and  it  is  noted  that  k0  =  k .  The  solution  given  in  equation  (4)  is  now 
converted  from  a  series  solution  into  a  closed-form  solution  using  an  analytical  expression  for 
the  summation  in  the  denominator.  Using  the  analytical  expression  given  by  (reference  1 1 ) 
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Equation  (5)  can  be  inserted  into  equation  (4),  and  this  yields  the  system  response 
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For  free-wave  propagation  to  exist  in  this  system,  the  denominator  of  equation  (6)  must  be  zero, 
i.e.. 
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The  term  k2  -  k2  is  not  a  frce-wavc  propagation  condition  for  the  reinforced  system  because 
when  k 2  =  k2  this  term  will  cancel  with  one  of  the  terms  in  the  numerator  of  equation  (6).  The 
relationship  between  wavenumber  and  free  wavenumber  (and  thus  frequency)  for  free-wave 
propagation  for  the  reinforced  system  is 
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and  this  yields  the  multi-valued  solution  to  the  wavenumber,  written  as 
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provided  that 
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which  can  be  alternatively  stated  as 
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Another  method  to  solve  this  problem  can  be  constructed  based  on  the  spatial  periodicity. 
Equation  (2)  can  be  changed  using  k„  =  k  +  2tt  n  /  L  so  that  an  infinite  number  of  indexed 
equations  are  written  as 
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These  individual  equations  can  be  placed  into  matrix  format  as 
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where  [D]  is  a  diagonal  matrix  whose  entries  are  equal  to  (A2  -  A,2 ) ,  [1]  is  a  matrix  whose  entries 
are  all  1,  {\v}  is  a  vector  of  the  unknown  displacements,  and  {f}  is  a  vector  of  all  ones  multiplied 
by  -  F0/T  .  The  matrix  in  equation  (12)  can  be  truncated  to  a  finite  number  of  terms  and  solved, 
yielding 
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Note  that  onee  equation  ( 1 3)  is  solved,  the  displacement  term  that  corresponds  to  the  system’s 
displacement  in  the  wavenumber  domain  (equation  (6))  is  w0(k0 ) .  If  enough  terms  are  chosen  in 
equation  (13),  this  solution  will  converge  on  the  expression  in  equation  (6).  For  the  system 
presented  in  equation  (12)  to  support  free  wave  propagation,  it  is  necessary  that 
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Therefore,  the  relationship  between  the  determinant  in  equation  (14)  and  the  system’s  free-wave 
propagation  requires  that 
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This  relationship  is  illustrated  graphically  in  figure  1.  For  this  example,  the  frequency  is 
275  Hz,  the  tension  is  100  N,  the  mass  per  unit  length  is  1  kg  m  ,  the  stiffener  spring  constant  is 
30,000  N/m,  the  stiffener  spacing  is  0.05  m,  and  37  terms  were  used  to  construct  the  determinant. 
The  solid  line  in  this  figure  is  a  plot  of  the  determinant  of  { [D]  +  [A.'  /(7X)][1] }  in  the  decibel 
seale  versus  wavenumber  from  -600  rad/m  to  600  rad/m.  The  solid  black  dots  are  the  free-wave 
propagation  wavenumbers  calculated  from  equation  (9),  and  these  match  the  minimum  values  of 
the  determinant  with  respect  to  the  wavenumber. 
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Figure  1.  Magnitude  of  Determinant  vs  Wavenumber  for  a  Second-Order  System  (solid  line) 
and  Free-Wave  Propagation  Wavenumbers  (dots)  Calculated  from  Equation  (9) 


3.  FOURTH-ORDER  SYSTEM 


The  motion  of  a  thin,  reinforced  plate  is  a  fourth-order  system  that  is  governed  by  the 
Euler-Bemoulli  plate  equation.  If  the  reinforcing  ribs  are  equally  spaced,  and  if  the  plate  has 
infinite  spatial  extent,  then  the  partial  differential  equation  modeling  of  this  motion  in  the  spatial¬ 
time  domain  is 
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where  ufx,/)  is  the  transverse  displacement  in  they-direction,  x  is  the  spatial  location  on  the 

plate,  t  is  time,  p  is  the  density,  h  is  the  height  of  the  plate,  f{x,t)  is  the  external  load  on  the 

plate,  K  is  the  stiffness  of  each  rib  per  unit  length,  L  is  the  distance  between  adjacent  ribs,  and  5 
is  the  Dirac  delta  function.  The  flexural  rigidity  of  the  plate  D  is 


D  = 


12(1  -u2)  ’ 


(18) 


6 


where  u  is  Poisson’s  ratio  and  E  is  Young’s  modulus  of  the  plate.  Using  a  line  load  that  is 
harmonic  in  time,  and  transferring  equation  (17)  into  the  wavenumber- frequency  domain,  yields 
(reference  10) 
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where  w(k)  is  the  transverse  displacement  in  the  v-direction,  k  is  the  wavenumber,  F0  is  the 
magnitude  of  the  applied  line  force  at  x  =  0,  and  Ay  is  the  free  (flexural)  wavenumber  of  the 
unreinforced  plate  and  is  given  by 


A,  = 


wv'4 

D  j 


(20) 


where  <y  is  the  frequency.  The  solution  to  this  equation  has  been  previously  solved  and  is  given  by 
vv(A)  - 1 
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The  solution  given  in  equation  (21)  is  now  converted  from  a  series  solution  into  a  closed  form 
solution  using  an  analytical  expression  for  the  summation  in  the  denominator.  Using  the 
analytical  expression  given  by  reference  1 1 


n~  cc 

Lis\n(kfL) 


cos 


(kL 


-  cosh 


V  2  J 


+  sinh(AAA) 


2  [kfL\ 

cos  — —  -cos 


kL 

2 


8A 


J  k,L  )  JkL' 
cos"  -  -cosy  — 


cos' 


J  kL 


-  cosh 


2  A/ 


J 


k,L A 


(22) 


Equation  (22)  can  be  inserted  into  equation  (21),  and  this  yields  the  system  response: 
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For  free-wave  propagation  to  exist  in  this  system,  the  denominator  of  equation  (23)  must  be  zero, 
i.e., 

A  =  0  ,  (26) 

where  A  is  rewritten  as 
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It  is  noted  that  Mace  (reference  8)  developed  a  similar  expression  for  this  relationship,  although 
his  parameter  set  was  non-dimensional.  The  term  k 4  -k4  is  not  a  free-wave  propagation 

condition  for  the  reinforced  system  because  when  k4  =k)  this  term  will  cancel  with  one  of  the 

terms  in  the  numerator  in  equation  (23).  Equation  (27)  yields  a  system  with  free-wave 
propagation  whenever 
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provided  that 
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Numerical  simulations  suggest  that  the  second  root  that  contains  a  minus  sign  inside  the  radieal 
is  extraneous.  If  the  conditions  listed  in  equation  (32)  are  satisfied,  then  the  wavenumber  where 
the  free-wave  propagation  occurs  can  be  calculated  using  the  formula 


k  =  ± 


-^-arceos(^)  + 


2n7t 

~r 


n  - 2, -1,0, 1,2  -  . 


(33) 


If  the  condition  listed  in  equation  (32)  is  not  satisfied,  then  the  system  does  not  support  free- 
wave  propagation  at  that  speeifie  frequency. 

An  alternative  method  to  solve  this  problem  ean  be  eonstrueted  based  on  the  spatial 
periodieity  using  equation  ( 1 9)  and  the  relationship  k„  =  k  +  2mi /  L  .  This  yields 


(k\ -kj)w.  i  (i'i)  =  "T"m^M'"^) 

( Dk A  -k}) vv0  ( k o )  =  ~ “■  "X  w>,  (k„)  •  (34) 

L )  L  n  T  , 

(Dk?  -  kAf)\V\  (k] )  =  -J7---pT  X  w"(k" ) 

L)  IJL  rj. 


The  individual  equations  from  equation  (34)  ean  be  placed  into  matrix  form  as 


[D]  + 


K 

DL 


m 


{w}  =  {f }  , 


(35) 


where  [D]  is  a  diagonal  matrix  whose  entries  are  (k*  -k*),  [1]  is  a  matrix  whose  entries  are  all 
1 ,  {w}  is  a  vector  of  the  unknown  displacements,  and  {fj  is  a  veetor  of  all  ones  multiplied  by 
-  Fn  /  D .  The  matrix  equation  in  equation  (35)  ean  be  truneated  to  a  finite  number  of  terms  and 
solved,  yielding 


w 


(  K  } 

D  +  — 

[1] 

{DL) 

(36) 
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Note  that  once  equation  (36)  is  solved,  the  displacement  term  that  corresponds  to  the 
displacement  of  the  system  in  the  wavenumber  domain  is  wa(k0).  If  enough  terms  are  chosen, 
this  solution  will  approach  the  expression  in  equation  (23).  For  the  system  presented  in  equation 
(36)  to  support  free-wave  propagation,  it  is  necessary  that 


det<[D]  + 


DL 


mho . 


(37) 


Therefore,  the  relationship  between  the  determinant  in  equation  (37)  and  free-wave  propagation 
of  the  system  requires  that 


det 


k\-k}  + 
K 


K 

DL 


DL 

K 

DL 


K 

DL 

ko  ~  kAf  + 

K 

DL 


K_} 
DL  j 


K 

DL 

K 

DL 

k\  —k}  + 


DL j 


=  0 


(38) 


if,  and  only  if, 

A' jsin(AyL) 
+  8  Dk) 


cos' 


rkL^ 


-cosh 


z  J 


k,L 


+  sinh(A'^L) 


COS' 


kfL  ^ 


cos2 

(kfL' 1 

1 

(  kL ) 

*> 

kL) 

-  cosh’ 

f  kfL\ 

-  cos~ 

COS“ 

l  2  J 

2  J 

2  J 

2  J_ 

v  2  j 
=  0 


-COS' 


kL^ 


*  J 


(39) 


This  relationship  is  illustrated  graphically  in  figure  2.  For  this  example,  the  following 
applies: 


•  Frequency  is  100  Hz. 

•  Thickness  of  the  plate  is  0.005  m. 

•  Young’s  modulus  of  the  plate  is  200e9  N  m  . 

•  Poisson’s  ratio  of  the  plate  is  0.3. 

•  Density  of  the  plate  is  7860  kg  m  \ 

•  Stiffener  spring  constant  per  unit  length  is  le6  N  m  . 

•  Stiffener  spacing  is  1.0  m. 

•  Eleven  terms  were  used  to  construct  the  determinant. 

The  solid  line  in  this  figure  is  a  plot  of  the  determinant  of  [D]  +  [K  /(DL)][  1]  in  the 
decibel  scale  versus  wavenumber  from  -20  rad/m  to  20  rad/m.  The  solid  black  dots  are  the  free- 
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wave  propagation  wavenumbers  calculated  from  equation  (33),  and  these  match  the  minimum 
values  of  the  determinant  with  respect  to  the  wavenumber. 


Figure  2.  Magnitude  of  Determinant  vs  Wavenumber  for  a  Fourth-Order  System  (solid  line) 
and  Free-Wave  Propagation  Wavenumbers  (dots)  Calculated  from  Equation  (33) 


4.  SUMMARY 


An  analytical  expression  for  free-wave  propagation  for  periodic  systems  has  been 
developed  by  solving  for  the  dynamic  behavior  using  two  separate  methods  and  equating  the 
results.  This  process  yields  a  relationship  between  a  diagonally-indexed,  full  determinant  of  a 
large  matrix  when  it  is  zero  and  a  simplified  analytical  expression.  This  process  was  undertaken 
for  both  a  second  and  fourth  order  system.  Two  example  problems  are  included  to  illustrate  the 
accuracy  of  this  process. 
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